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lLX^vo},  abstract 


This  thesis  considers  the  influence  of  four  secondary  effects 
on  the  analysis  of  a  simple  ell-shaped  piping  structure  when  it  is 
uniformly  heated.   The  four  effects  are:   (a)   axial  deformation, 
(b)   shearing  deformation,   (c)   beam-column  effect,  and   (d)   dif- 
ference between  arc  and  chord.   They  are  usually  neglected  in  a 
piping  flexibility  analysis. 

In  a  paper  of  several  years  ago,  J.  E.  Brock  developed  a  theory 
for  this  investigation  but  diA|not;:$^jg^it  any  numerical  results  or 
conclusions.   The  present  thesis  reviews  Brock's  theory  and  corrects 
some  errors  in  coefficients.   It  then  describes  a  digital  computer 
program  for  IBM  System  360  which  performs  the  corresponding  calculations. 
Finally,  it  draws  the  new  and  quite  significant  conclusion  that  con- 
ventional piping  stress  analysis,  which  neglects  these  secondary  effects, 
can  result  in  gross  errors  in  evaluating  stresses  in  piping  config- 
urations which  are  likely  to  be  encountered  in  piping  practice. 
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SECTION  I 
INTRODUCTION 

1-1  Background 

In  1964  Professor  John  E.  Brock  of  the  Mechanical  Engineering 
Department  of  the  Naval  Postgraduate  School,  Monterey,  California, 
published  a  paper  entitled,  "Some  Secondary  Effects  in  a  Simple 
Piping  Structure  Under  Heating". 

The  secondary  effects  referred  to  are  axial  deformation,  shear, 
beam  column  action,  and  difference  between  arc  and  chord.   The 
structural  element  considered  is  an  ell-shaped  piping  configuration 
with  a  square  corner. 

In  the  normal  analysis  of  such  a  structure,  the  above  mentioned 
secondary  effects  are  neglected,  being  regarded  as  insignificant  in 
relieving  the  stresses  caused  by  thermal  expansion.   Only  the 
principal  sources  of  relief,  bending  and  torsion,  are  usually  consid- 
ered. 

Professor  Brock  undertook  to  make  a  complete  flexibility  analysis 
of  this  piping  configuration  including  all  of  the  four  secondary 
effects.   The  1964  paper  reports  his  results.   It  also  reports  that 
a  computer  program  had  been  written  to  investigate  this  analysis 
and  the  resulting  system  of  equations.   The  program  had  been  written 
but  was  not  developed  sufficiently  to  thoroughly  investigate  the 
relationships  involved. 


John  E.  Brock,  "Some  Secondary  Effects  In  A  Simple  Piping  Structure 
Under  Heating,"  Journal  of  Applied  Mechanics,  Vol  31,  March  1964, 
pp.  88-90. 


1-2  Objectives  of  Work 

The  object  of  this  thesis  was  to  write  a  computer  program  for 
a  large  digital  computer  that  would  be  capable  of  performing  the 
calculations  required  by  Professor  Brock's  theoretical  development 
and  one  that  would  allow  application  of  the  theory  to  any  engineering 
problem  that  involves  a  piping  structure  of  the  given  configuration. 
Further,  it  attempts  to  determine  under  what  conditions  these 
secondary  effects  may  become  significant  in  calculation. 
1-3  Assumptions  and  Limitations 

The  investigations  that  are  reported  herein  are  based  on  an  ell- 
shaped  piping  configuration  with  a  square  corner.   No  other  shape  is 
considered.   Also,  the  configuration  is  assumed  to  have  deformation 
only  in  the  plane  of  the  legs  a  and  b. 

It  is  further  assumed  that  deformations  and  displacements  are 
small  compared  to  the  longitudinal  dimensions  of  the  pipe,  except 
that  a  first  approximation  is  included  for  the  difference  between 
arc  length  and  chord  length  of  the  deformed  pipes. 
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1-4  Notation 

TABLE  I 

NOMENCLATURE 

a,b  length  of  straight  pipe  elements 

A  cross  sectional  area  of  pipe  material 

A  cross  sectional  area  of  pipe  contents 
r 

e  unit  thermal  strain 

E  Young's  modulus  of  elasticity 

F. ,F  structural  forces  in  ell-configuration 

g.  .  coefficients  defined  in  equations  (16)-(19) 

G  shearing  modulus  of  elasticity 

H  (l-ZJJ 0/E  (l-  pz)  =  pressure  coefficient 

I  moment  of  inertia  of  pipe  cross  section 

L  length  of  pipe 

M,N  bending  moments 

p  internal  pressure 

P  lateral  force 

Q  axial  compressive  force 

u  ur* 

v  (-us)'"* 

VJ  QL2/EI 

x  axial  coordinate 

X  horizontal  deflection  of  joint 

y  lateral  coordinate 

Y  vertical  deflection  of  joint 

°\  coefficients  (see  Table  II) 
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TABLE  I  (cont.) 

Yl   .  P^  coefficients  defined  in  equations  (15A) 

and  (15B) 
0  yf/Afj  =   effective  shearing  strain 

tip  axial  extension 
A  tip  lateral  deflection 

G>ij   £-2  conveniently  small  positive  numbers 

S  shear  distribution  factor  (see  equation  (9)) 

T)  ^EI/AGL2 

C?  rotation  of  joint 

A  difference  between  arc  and  chord  (first 

approximation) 
Ml)  Mi.  coefficients  in  equations  (9),  (10),  and  (11) 

Poisson's  ratio 

(inside  diameter  of  pipe) /(outside  diameter 
of  pipe) 

tip  rotation 


r 


tj  xu/L 
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SECTION  II 
THEORETICAL  ANALYSIS  AND  TECHNIQUE  FOR  INVESTIGATION 


2-1  Review  of  Flexibility  Analysis  for  the  Ell-Shaped  Pipe 

Before  presenting  the  investigation  of  this  ell-shaped  piping 

configuration   in  detail,  it  is  felt  that  certain  principles  re- 
garding the  analysis  merit  some  review. 
If  the  tip  loaded 

cantilever  beam  shown  in 

Figure  1  is  analyzed  by 

writing  a  moment  equation 

about  point  A  the  following 

equation  results: 


Figure  1 


(1) 


Ely"  =  m  +  Q(a-y)+  p(l--x) 


Solving  this  differential  equation  and  applying  boundary 
conditions  yields 
(2) 

EIu? 


s   yields 


for  boundary  conditions 

a. 

y(0)  =  o 

b. 

yd)  =  A 

c. 

y'(0)  =  V 

d. 

y'(L)  =  Cf 

See  Table   I    for  nomenclature. 


In  addition 


(3) 


A  =  LMd,  +  PL(o(a  +  1Jof3)](LVEIN) 
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and 


W      cf>  »  [Mc*j+  PL(o(>+T)c<4)](L/El) 


The  alpha  coefficients  appearing  in  the  above  and  subsequent 

equations  are  defined  in  Table  II  with  alternate  forms  to  facilitate 
computation. 

In  terms  of  all  four  secondary  effects, &  is  given  by  equation  (5) 

(5)  5  =  el  -  afoL/  _    (Q-A#)l  _    7\ 

The  first  term  on  the  right  hand  side  of  equation  (5)  represents 
thermal  expansion,  the  second  represents  axial  contraction  resulting 
from  hoop  tension  caused  by  internal  pressure,  the  third  represents 
axial  contraction  due  to  axial  compression  in  the  metal,  and  the 
last  represents  the  effective  axial  contraction  resulting  from  the 
difference  between  arc  and  chord.   The  pressure  terms  may  be  combined 
to  yield  the  equation, 

(6)  £  =    (e+f>n)l_  -    Qi_/ae  -   7\ 

where  A,   is  given  by  , 

when  higher  order  terms  are  neglected  (see  Appendix  G) .   If  alpha 
coefficients  are  substituted, 

(8)    /\  =  ( M'o(8  +  2.0  MPLdj  +  P'L^X  I? A??) 
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factor 

(9) 


Other  pertinent  relationships  include  the  shear  distribution 
2 


where 


(10) 


and 


(1 


/*,  *  (14-  i4-i)+8J)/<o(/+tjy 


d  ^ =  (io+zoi>t 5v»9/3 o i-py 


To  apply  the  foregoing  analysis  to  an  ell-shaped  pipe,  the 
corresponding  quantities  are  defined  as  they  appear  in  that 
configuration.   The  basic  structure  is  that  shown  in  Figure  2.   If 
the  configuration  is  subjected  to  uniform  heating,  thermal  expansion 
creates  a  situation  resembling  that  of  Figure  3. 


/////. 


T 


///// 


h* *-t 

Figure   2 

The  resulting  free  body  diagrams  for  the  two  legs  a  and  b  is  depicted 
in  Figure  4. 


John  E.  Brock,  "Shear  Distribution  in  Piping,"  Heating  Piping  and 
Air  Conditioning,  Vol  35,  January  1963,  pp.  141-143. 
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ueu/ 


Figure  4 


Referring  to  the  horizontal  leg  as  part  1  and  to  the  vertical  leg  as 
part  2,  Table  III  list  the  quantities  corresponding  to  the  general 
case  for  each  part. 


TABLE  III 

ITEM 

PART  1 

PART  2 

L 

a 

b 

P 

-Fi 

F2 

Q 

F2 

Fl 

M 

N 

-N 

A 

-Y 

X 

6 

X 

Y 

f 

e 

9 

Geometrical  continuity  requires  equating  corresponding  expressions 
for  X,Y,  and  9;  i.e., 


6k  =    Ab 
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or 


(12)  ea -  ^  -  ft  +  M  = [-4a,  +  F4b (av+  Sjefe'j)  b/fci 


A  superior  bar  indicates  that  the  quantity  has  been  evaluated 
using  the  parameters  of  part  2  and  no  superior  bar  indicates  evaluations 
with  the  parameters  of  part  1. 

Solving  equation  (14)  for  the  moment,  N,  yields 

(15) 

where    (15A) 


and    (15B) 


Now   define; 
(16) 

(17) 

(18) 
(19) 


and  the  system  of  equations  become 


(20) 


Sin   3.i 

<3*l  <i»i 


F, 


(e+pw)a  -  A 
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In  order  to  solve  for  the  inknowns  F  and  F  ,  an  iterative 
procedure  is  strongly  suggested  by  the  nature  of  the  equations. 

By  assuming  initial  values  of  F1  and  F_, which  can  be  any 
value,  one  can  calculate  finite  values  for  all  the  coefficients 
g.  .  and  the  quantities  on  the  right  of  equation  (20)  .  This  equation 
is  solved  for  new  values  of  F,  and  F„  which  in  turn  provide  improved 
values  for  the  next  iteration.   This  procedure  of  iteration  is 
continued  until  satisfactory  convergence  is  reached.   Once  convergence 
is  obtained  such  quantities  as  moment,  0,  and  other  desired  values 
are  obtained  using  the  converged  values  of  F..  and  F  . 
2-2  Dimensionless  Parameters 

The  large  number  of  variables  involved  in  this  problem  could 
produce  millions  of  different  problems.   In  order  to  pursue  the 
investigation  as  thoroughly  as  possible  and  at  the  same  time  retain 
some  compactness  in  the  results,  four  dimensionless  parameters  were 
used.   These  are: 

(1)  Leg  Ratio  a/b 

(2)  (Radius  of  Gyration) /(Total  Length) 

(3)  Pressure/Modulus  of  Elasticity 

(4)  Thermal  Strain 

To  cover  reasonable  ranges  of  all  four  parameters  evidently 
requires  solution  of  a  large  number  of  problems  each  of  which  is 
itself  rather  demanding.  Accordingly,  use  of  a  large  digital  computer 
is  definitely  required. 

Each  of  the  four  secondary  effects  can  either  be  included  or 

4 
excluded  in  the  calculation.  This  produces  2   or  16  different  possible 

combinations  of  secondary  effects. 
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2-3  Requirements  Imposed  on  Digital  Computer  Program 

A  computer  program  for  studying  this  problem  has  certain  basic 
requirements . 

A.  First  of  all, it  must  be  capable  of  accepting  certain  input 
data,  which  includes  a  basic  physical  description  of  the 
problem  as  well  as  control  variables  that  instruct  the 
program  in  how  to  treat  these  data. 

B.  From  the  input  information  it  must  recognize  the  parameter 
to  be  investigated  and  must  load  the  various  values  of  the 
parameter  over  the  desired  range. 

C.  For  each  parameter  value  in  the  given  range  it  must: 

1.  Solve  for  the  F's  by  the  iteration  procedure  and  then 
determine  other  desired  values. 

2.  Accomplish  C.l.  for  all  different  combinations  of  the 
four  secondary  effects. 

D.  It  must  increment  the  parameter  value  and  repeat  step  C 
until  results  for  all  values  in  the  range  have  been  found. 

E.  It  must  provide  results  in  either  tabular  or  graphical  form, 
The  following  section  describes  the  computer  program  in  more 

detail  and  shows  how  the  above  listed  demands  are  met. 
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SECTION  III 
DESCRIPTION  OF  COMPUTER  PROGRAM 

3-1  General  Remarks 

As  mentioned  in  Section  II  and  suggested  in  reference  2,  a 
computational  procedure  is  evident.   Assume  initial  values  of 
F1  and  F  (which  may  be  taken  to  be  zero),  evaluate  the  coefficients 
g. .  and  the  terms  A  and  A  and  upon  salving  the  system  of  equations  (20), 
obtain  improved  values  of  F1  and  F_,  continuing  until  satisfactory 
convergence  is  reached.   This  idea  along  with  the  other  basic  require- 
ments (A-E  of  Section  II)  composes  the  foundation  upon  which  the 
computer  program  is  constructed. 

The  program  was  first  written  in  Fortran  60  computer  language 
for  use  on  the  Control  Data  Corporation  Model  1604  digital  computer 
located  at  the  Naval  Postgraduate  School,  Monterey,  California.   In 
the  midst  of  the  investigation  reported  in  this  paper,  the  CDC 
computer  was  replaced  at  the  Postgraduate  School  by  an  International 
Business  Machine  Model  360/67.   The  present  version  of  the  program 
is  written  in  Fortran  IV  language  for  use  on  the  latter  machine. 

All  computations  in  the  program  are  executed  using  the  Double 
Precision  Capability  of  the  system.   This  feature  improves  the 
accuracy  of  the  results  but  unfortunately  requires  more  time  than 
the  single  precision  calculation. 

An  effort  was  made  to  keep  the  program  as  general  and  as  flexible 
as  possible.   All  of  the  information  needed  to  specify  a  given  problem 
is  read  directly  from  data  cards.   The  information  on  the  data  cards 
defines  the  physical  condition  to  be  investigated  and  gives  the  program 
additional  control  information  concerning  the  parameter  to  be  used 
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in  evaluation,  range  of  the  parameter,  combinations  of  the  secondary 
effects  to  be  investigated,  convergence  criteria,  etc. 

It  allows  the  user  four  different  parameters  with  which  to 
evaluate  the  problem.   It  allows  for  all  sixteen  possible  combinations 
of  the  four  secondary  effects  and  computes  all  sixteen  of  these  for 
each  value  of  the  parameter. 

It  has  optional  output.   In  the  data  cards  the  user  specifies 
the  desired  type  of  output:   tabular,  graphic,  or  both. 

To  describe  the  computer  program  more  completely  requires 
explaining  in  more  detail  the  functions  of  its  seven  major  parts. 

A.  Input 

B.  Loading  of  Parameter  Values 

C.  Consideration  of  All  16  Different  Combinations 

D.  Computation  by  Iteration 

E.  Incrementing  Parameter  to  Next  Value  in  Range 

F.  Output 

G.  Subroutines 
3-2   Input 

The  term  input  is  used  here  in  referring  to  the  information  that 
is  supplied  to  the  program  via  the  data  cards.   In  all,  there  are 
fifty-five  different  items  read  in  as  data.   These  items  completely 
define  the  physical  situation  under  consideration.   Further,  it  informs 
the  program  of  the  dimensionless  parameter  to  use  in  the  investigation, 
the  range  of  the  parameter,  conditions  to  investigate,  and  so  on. 
Each  of  the  variables  read  in  as  data  is  listed  below  along  with  a 
description  of  its  function  in  the  program.   The  order  of  this  listing 
is  the  order  in  which  they  appear  on  the  data  cards.   The  names  of  the 
variables  are  the  names  which  are  used  as  input.   The  same  physical 
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quantity  may  be  called  by  other  names  later  in  the  program  and  it 
should  be  remembered  that  the  variable  name  mentioned  below  is 
applicable  only  in  the  input  section  of  the  program. 


R 

S 

ZIN 

ZOUT 

DTEMP 

Q 
E 
COEFI 

KPARl 


length  of  the  shorter  leg,  a,  (inches) 
length  of  the  longer  leg,  b,  (inches) 
pipe  inside  diameter  (inches) 
pipe  outside  diameter  (inches) 

temperature  difference  (degrees  fahrenheit)  = 

(T      -  T       ) 
final    initial 

internal  pressure  (psi) 

modulus  of  elasticity  of  pipe  material  (psi) 

coefficient  of  linear  expansion  of  pipe  material 
(in/in  -  °F) 

Integer  variable  which  can  assume  the  integer  values 
1  to  4  inclusive,  that  designates  the  parameter 
to  be  used  in  the  investigation  as  follows: 

KPARl  =  1  =  Leg  Ratio 

KPARl  =  2  =  (Radius  of  Gyration) /(Total  Length) 

KPARl  =  3  =  Pressure/Modulus  of  Elasticity 

KPARl  -  4  =  Unit  Thermal  Strain 

Integer  variable  that  specifies  the  size  of  the  major 
parameter  increment  (see  Section  3-3).  This  can  be 
any  value  but  for  graphic  output  should  not  require 
more  than  thirty  stored  values  of  the  parameter. 


23 


NN  Integer  variable  that  designates  the  initial  value 

of  N,  itself  an  integer  variable  which  is  used 
as  a  subscript  for  the  parameter  values.   This 
(the  value  of  NN)  is  useful  if  the  program  ever 
exceeds  the  time  allotted  for  one  run  and  it  is 
terminated  by  the  operator,  or  if  the  given  size 
of  increments  do  not  yield  convergence  over  the 
desired  range.   By  giving  NN  the  value  of  the  last 
value  of  N  on  the  terminated  run,  the  user  is  able 
to  modify  certain  other  control  variables  such  as 
L,  LL1,  FTE's,  etc.,  for  the  next  computer  run  and 
begin  computation  at  the  point  where  the  last  run 
ended,  regardless  of  the  reason  for  termination. 
For  the  first  run  of  a  given  parameter  or  problem 
the  value  of  NN  is  always  1. 

MA  Integer  variable  that  designates  the  final  value  of  the 

subscript  N  in  the  parameter  range. 

MB  Integer  variable  which  designates  the  initial  value  of 

MT .   The  significance  of  MT  is  explained  in  more 
detail  in  the  section  defining  loading.   MB  can  take 
on  any  value  from  1  to  16. 

MC  Integer  variable  that  indicates  the  final  value  of  MT. 

MC  >MB  always  and  1<MC<£16.   Both  of  the  above 
conditions  must  be  satisfied.   The  program  is  designed 
in  such  a  way  that  if  MT  does  not  equal  its  final 
value,  MC,  then  MT  will  be  incremented  by  1 ,  or  the 
next  combination  of  secondary  effects  will  be 
calculated  using  the  same  parameter  value.   What 
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this  in  effect  does  is  have  the  computer  make  all  the 
calculations  for  all  the  combinations  between  and 
including  MT  =  MB  and  MT  =  MC .   For  example,  if 
MB  =  3  and  MC  =  7  for  a  given  set  of  input  values, 
MT  would  take  on  values  3,  4,  5,  6,  and  7  for  each 
value  of  the  parameter.   Consequently,  the  output 
would  reveal  results  for  the  combination  correspond- 
ing to  the  5  different  values  of  MT. 

MD  Integer  variable  that  designates  the  type  of  output 

desired . 
1   MD   3 

MD  =  1  Tabular  output  only 
MD  =  2  Graphic  output  only 
MD  =  3    Both  tabular  and  graphic  output 

MIT  Integer  variable  that  indicates  the  limit  to  the 

number  of  iterations  the  program  will  perform.   If 
this  number  is  reached  before  the  convergence 
criterion  is  met,  the  program  will  terminate 
operation  on  its  current  set  of  data,  give  the 
necessary  output,  then  proceed  to  the  next  set 
of  data. 

KKK  Integer  variable  that  assigns  a  value  to  the  constant 

IA.  IA  is  used  in  a  relationship  that  determines 
the  improved  values  of  F's  to  be  used  in  the  next 
iteration  on  the  basis  of  the  last  assumed  values 
of  the  F's.   The  relationship  involved  is: 

4  +  IA 
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LL1 


EP2 


EXSCAL 


YSCALE 
IXUP 


IYRIGH 


FTE11-FTE162 


and  a  similar  relationship  exists  for  the  improved 
F2.   It  was  found  that  generally  as  the  value  of  N 
increases  IA  should  increase  also,  thereby  making 
the  improved  value  closer  to  the  old  value  than 
the  new  value. 

Integer  variable  that  specifies  the  size  of  the  sub- 
increments  of  the  parameter.   (See  section  concerning 
loading) . 

This  value  sets  the  convergence  criteria.  After  each 
iteration  the  new  values  of  the  F's  are  compared  to 
the  old  ones  using  the  equation 

V  (F1(WEW))*-+-  (F^N/ew)" 

If  this  is  satisfied  the  values  are  considered  converged, 

This  quantity  is  the  argument  used  in  the  DRAW  subroutine 
(see  Appendix  A-4)  that  assigns  the  scale  to  the  X-axis 
or  abscissa  axis  when  graphic  output  is  desired. 

Same  a  EXSCAL  above  except  it  governs  scale  on  Y-axis. 

Another  argument  for  the  DRAW  subroutine  that  specifies 
how  many  inches  up  the  graph  one  wants  the  X-axis 
located. 

Specifies  how  many  inches  from  the  left  side  of  the  graph 
one  wants  the  Y-axis. 

The  32  different  values  are  the  last  converged  values 
of  Fl  and  F2  for  a  particular  combination  of  secondary 
effects.   The  final  digit  distinguishes  between  Fl 
and  F2 .   The  other  numerical  digits  run  from  1  through 
16  and  correspond  to  values  of  MT.   These  FTE  values 
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are  retained  by  the  computer  and  used  as  the  first 
approximation  for  the  F's  in  the  corresponding 
combination  when  the  parameter  has  been  incremented. 
In  most  investigations  the  values  assigned  by  the 
data  card  are  all  zero.   These  values,  however,  are 
just  the  initial  values  assigned  to  the  FTE's.   Once 
computation  begins  the  computer  takes  over  the  task 
of  assigning  these  values.   Every  time  convergence  is 
reached  the  values  of  two  FTE's  change  and  as 
convergence  is  reached  for  all  combinations  and  one 
parameter  value,  all  32  change  their  values.   Use  of  the 
variables  aids  the  computation  process  in  two  ways. 
The  nature  of  the  governing  equations  causes  difficulty 
with  convergence  under  some  conditions.   Therefore, 
if  the  first  assumed  values  of  the  F's  are  near  the 
correct  values  a  better  chance  for  convergence  exists. 
If,  on  the  other  hand,  the  guess  is  arbitrary,  or  not 
close  to  the  converged  value,  the  system  of  equations 
may  not  converge.   If  the  last  converged  values  are 
assumed  as  the  first  approximations  with  the  incremented 
parameter,  they  should  be  good  estimates  depending  on 
the  size  of  the  parameter  increment.   The  smaller  the 
parameter  increment  the  better  the  estimate  should  be. 

Another  utilization  of  the  32  variables  allows 
the  user  the  capability  of  picking  up  computation  at 
any  point  in  a  parameter  range.   Suppose,  due  to  excess 
time  the  run  was  terminated  before  the  entire  range 
was  completed.   Then,  on  the  subsequent  run,  the  user 
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would  be  able  to  start  at  the  same  place  in  the  parameter 
range  that  the  previous  run  ended.   This  provides 
continuity  without  duplication  and  the  ability  to  read 
reasonable  first  approximations  to  the  F's  in  as  data 
regardless  of  where  in  the  range  one  desires  to  continue. 

3-3  Loading  of  Parameter  Values 

The  integer  variable  KPAR1,  read  from  the  data  card,  instructs 
the  program  which  dimensionless  parameter  will  be  used  to  investigate 
the  structure.   For  a  given  set  of  data  (one  set  of  six  data  cards) 
only  one  parameter  can  be  investigated.   Each  of  the  four  values  that 
KPAR1  may  take  on  indicates  a  different  loading  scheme  for  the  six 
subscripted  variables.   Once  the  program  knows  the  value  of  KPAR1  it 
proceeds  to  the  proper  loading  scheme  for  that  parameter. 
The  six  subscripted  variables  involved  are: 

AA(N)  leg  a 

BB(N)  leg  b 

XDIN(N)        pipe  inside  diameter 

XDOUT(N)       pipe  outside  diameter 

XP(N)  internal  pressure 

ALIT(N)        unit  thermal  strain 
Here  again  the  variable  names  associated  with  the  physical 
quantities  apply  only  in  the  loading  scheme.   The  same  physical  quantity 
may  have  another  variable  name  at  another  place  in  the  program.   See 
Appendix  A  for  a  complete  listing  of  computer  program  nomenclature. 

The  subscript  N  takes  on  integer  values  from  1  to  MA.   For  each 
value  of  N,  the  six  subscripted  variables  have  values  stored  in  the 
computer  memory.   In  each  of  the  different  loading  schemes  only  one 
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of  the  six  subscripted  variables  takes  on  different  values  for  each 
different  value  of  N.   The  other  five  have  fixed  values. 

Associated  with  N  in  the  loading  procedure  are  two  other  integer 
variables,  MA  and  L.   The  values  of  these  two  quantities  determine 
the  range  of  the  parameter  and  the  size  of  the  main  parameter  increment. 

An  important  feature  of  the  loading  process  is  that  the  parameter 
values  are  loaded  in  such  a  way  as  to  produce  the  least  critical  loading 
conditions  on  the  structure  for  the  parameter  value  associated  with  N  =  1  to 
the  most  critical  loading  condition  when  N  =  MA.   This  essentially, 
means  that  at  the  least  critical  end  of  the  range  the  initial  assumption 
of  the  arbitrary  values  of  the  F's  will  probably  yield  convergence.   How- 
ever, as  the  parameter  proceeds  to  the  more  critical  end  of  the  range, 
the  facility  utilizing  the  FTE's  becomes  more  and  more  necessary  in 
achieving  convergence. 
3-4  Combination  of  Secondary  Effects 

The  capability  that  allows  the  computer  program  to  solve  the  system 
of  equations  for  all  sixteen  different  combinations  of  secondary  effects 
for  each  parameter  value  is  afforded  through  the  subroutine  FACTOR. 
To  supply  the  subroutine  with  information  it  needs,  the  integer  variable 
MT  is  used.  MT  takes  on  values  1-16  for  each  parameter  value  in  the 
range.   Each  of  the  16  values  represents  a  different  combination  of 
secondary  effects.   When  the  subroutine  is  called  with  MT  specified, 
it  responds  by  supplying  the  main  program  with  four  factors  and  a 
particular  designator.   These  four  factors  are  multiplied  by  the  terms 
in  the  system  of  equations  that  represent  the  four  different  secondary 
effects.   Each  of  the  factors  can  have  only  two  different  values:   1.0 
and  0.0.   It  is  therefore  evident  that  if  the  combination  requires  the 
suppression  of  an  effect  in  the  calculation,  the  governing  factor  is 
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zero;  on  the  other  hand  if  the  effect  is  to  be  included  the 
governing  factor  is  unity. 

The  main  program  is  so  designed  that  all  of  the  sixteen  different 
combinations  do  not  have  to  be  calculated  each  time,  say  in  a  situation 
where  the  user  is  only  interested  in  one  or  two  combinations.   If  such 
is  the  case  the  integer  variables  MB  and  MC  are  manipulated  to  provide 
the  combinations  desired.   MB  sets  the  initial  value  of  MT.   After 
each  combination  is  completely  solved  for  all  desired  quantities  the 
main  program  tests  the  value  of  MT  against  that  of  MC.  MC  represents 
the  final  value  of  MT  to  be  considered  for  each  parameter  value  and 
in  the  test  if  MT  is  less  than  MC  its  value  is  incremented  by  1  and 
that  combination  is  investigated.   This  process  continues  until 
MT  =  MC  and  after  the  MT  =  MC  combination  is  solved  the  program  increments 
the  parameter.   The  limitation  on  this  capability  lies  in  the  fact  that 
combinations  represented  by  values  of  MT  between  MT  =  MB  and  MT  =  MC 
will  be  calculated  even  if  some  of  these  intermediate  cases  are  not 
specifically  wanted. 

The  designator  referred  to  earlier  is  realized  as  a  label  in  the 
printout  of  the  results.   It  is  printed  in  a  certain  column  of  the 
results  and  is  a  four  character  identifier  that  tells  the  value  of 
each  factor  in  that  particular  combination  (see  the  explanation  of  this 
subroutine  for  specific  meanings  for  these  designators  in  Section  3.8). 
3-5   Computation  by  Iteration 

The  iterative  procedure,  as  mentioned  earlier,  is  that  of  assuming 
initial  values  of  Fl  and  F2  (which  may  be  taken  to  be  zero),  evaluating 
the  coefficients  e. .  and  the  terms  /\  and  /}  and  upon  solving  the  system 
of  equations  (20),  obtaining  improved  values  for  the  F's.   This  is 
continued  until  satisfactory  convergence  is  reached.   Through  the  use 
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of  the  values  of  the  FTE's  an  attempt  is  made  to  minimize  the  number 
of  iterations  required.   This  saves  computer  time  as  well  as  increases 
the  probability  of  convergence. 

After  each  iteration  the  newly  calculated  values  of  the  F's  are 
tested  against  the  old  ones  using  the  following  convergence  criterion: 

|Fi(rttw)  -  F1(old)|+  |Fa(MEW)-  P2(M| 

If  this  criterion  is  not  satisfied,  the  program  selects  improved 
values  for  the  next  iteration  and  sends  these  newly  selected  assumptions 
of  Fl  and  F2  through  the  system  of  equations  again.   The  process  for 
the  selection  of  the  new  F's  is  defined  by  the  following  two  relation- 
ships : 

F1MFXT=  4*Flti  +  lA*Fl  .      FfcNEKT8   4»*F2N  +IA*  F£ 
4-+IA         >  4  +  XA 

where  FINEXT  and  F2NEXT  are  the  improved  values  of  Fl  and  F2  respectively 
to  be  used  in  the  next  iteration.   FIN  and  F2N  are  the  new  values  of 
F.  and  F   just  computed  and  Fl  and  F2  in  the  equations  are  the  values 
of  F.  and  F  assumed  at  the  beginning  of  the  iteration.   The  factor 
IA  is  an  integer  whose  value  is  set  by  the  integer  variable,  KKK, 
read  in  from  the  data  card.   It  is  obvious  that  the  improved  value  can 
be  closer  to  the  older  value  or  the  newer  value  at  the  discretion  of 
the  user.   The  nature  of  the  equations  is  such  that  when  the  parameter 
value  produces  less  critical  loading  conditions,  convergence  is  achieved 
more  rapidly  if  the  improved  value  is  taken  to  be  closer  the  new  values. 
Conversely,  when  the  parameter  is  in  a  more  critical  region,  a  selec- 
tion of  an  improved  value  closer  to  the  old  value  produces  convergence 
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more  rapidly.   The  term  "critical  region"  refers  to  that  range  of  the 
parameter  where  the  system  of  equations  behave  in  a  highly  nonlinear 
fashion. 

There  is  no  set  rule  for  this  selection  of  improved  values  and 
the  best  instruction  is  gained  through  experience  with  the  program. 
Using  values  of  IA  between  0  and  10  is  normally  adequate  for  most 
situations. 

If  the  convergence  criterion  is  satisfied,  the  latest  calculated 
values  of  the  F's  are  accepted  as  the  correct  ones  for  the  given  value 
of  the  parameter  and  the  given  combination  of  secondary  effects. 
Using  these  values  of  the  F's  other  necessary  quantities  are  calculated. 
The  program  then  goes  to  the  next  combination  of  secondary  effects 
and  begins  iteration  for  the  solution  of  the  F's  for  that  situation. 
The  first  approximations  for  the  F's  are  the  converged  values  found 
for  that  combination  with  the  previous  parameter  value.   When  the 
program  has  solved  the  equations  for  all  combinations  desired,  it 
increments  the  parameter  and  the  process  starts  all  over  again. 
3-6  Incrementing  the  Parameter 

Two  different  parameter  increments  are  utilized  in  the  program. 
They  are  referred  to  hereafter  as  &PAR  and  Qpar.   The  two  can  be 
described  as  follows: 

If  PAR(N)  is  the  parameter  for  a  particular  value  of  the  subscript, 
then  PAR(N  +  1)  =  PAR(N)  +  &PAR.   Therefore,  A  PAR  is  the  increment 
whose  size  is  the  difference  between  two  consecutive  stored  parameter 
values.   The  size  of  APAR  is  determined  by  appropriate  input  data 
and  by  the  control  variable  L  (see  description  of  L  in  Section  3-2). 
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The  sub- increment ,  <5par  is  used  to  give  the  program  the  capability 
of  making  the  parameter  as  small  as  deemed  necessary  in  order  to  achieve 
convergence.   The  basic  theme  of  the  investigation  of  any  of  the  four 
dimensionless  parameters  is  to  proceed  from  that  end  of  the  range 
that  provides  the  least  critical  load  toward  the  end  which  produces 
the  most  critical  condition  in  the  range.   Each  time  the  parameter 
is  incremented,  the  last  converged  values  of  the  F's  are  used  as 
the  first  approximation.   By  making  the  increments  smaller,  the  first 
approximation  is  closer  than  if  the  parameter  had  a  larger  increment. 
Therefore,  it  is  conceivable  that  by  making  these  increments  small 
enough,  convergence  can  always  be  realized.   It  is  for  this  reason 
that  the  program  uses  increments  of  two  different  sizes. 

For  the  graphic  output  capability,  it  is  necessary  to  store  many 
values  of  the  variables  that  are  used  in  the  graphs.   The  storage 
capacity  of  the  computer  thus  establishes  a  lower  limit  for  the  larger 
increment.   In  order  to  further  decrease  the  size  of  the  increment  the 
smaller  increment,  <SPa^»  is  used.   Except  when  they  coincide  with 
APAR  values,  parameter  values  resulting  from  this  smaller  increment 
and  the  corresponding  calculated  values  are  not  stored,  but  used  only 
to  obtain  new  FTE '  s  for  the  subsequent  parameter  value.   Another  way 
of  expressing  the  process  is  to  say  that  after  all  necessary  calculations 
are  made  for  one  parameter,  the  parameter  value  is  incremented  by  £par. 
Using  the  last  converged  values  of  the  F's  as  first  approximations, 
converged  values  for  the  new  parameter  value  are  sought.   Thus  the 
incrementation  continues  over  the  parameter  range.   If  the  first 
parameter  value  considered  was  a  stored  value,  then  after  LL1  sub- 
increments  the  next  stored  parameter  value  will  be  reached.   LL1  is 
the  integer  variable  read  from  the  data  card. 
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Figure  5  aids  in  understanding  the  use  of  the  two  increments. 

ill  Sahincrtroents 
^3 


^(P&rimeHr) 


PARAMETER. 

Figure  5 

From  Figure  5  the  following  relationships  are  apparent: 

LL1  x  £par  =  APAR        PAR(N  +  1)  =  PAR(N)  +  APAR 
Also,  for  the  parameter  values  corresponding  to  the  solid  vertical 
lines,  necessary  calculated  values  are  stored  for  graphic  or  tabular 
output.   For  parameter  values  corresponding  to  the  broken  vertical 
lines  only  converged  values  of  the  F's  are  calculated  for  use  as 
the  first  approximation  when  the  parameter  is  incremented. 

To  summarize,  each  time  the  parameter  is  incremented  it  is 
increased  by  the  amount  ^par.   Every  LL1  times,  all  calculated  values 
necessary  for  output  are  stored.   The  size  of APAR  is  somewhat 
restricted  by  the  storage  capacity  of  the  computer  but  the  size  of 
5  par  has  no  restrictions.   It  must  be  remembered,  however,  that  by 
decreasing  the  size  of  ©par,  tne  computing  time  for  a  given  range  of 
the  parameter  is  increased  accordingly.   Therefore,  in  the  interest 
of  time  on  the  computer,  the  value  of  £par  should  be  as  large  as  possible, 
commensurate  with  providing  the  necessary  convergence. 
3-7   Output 

The  user  has  some  choice  concerning  the  type  of  output.   The 
results  can  be  represented  in  tabular  form,  graphical  form,  or  both. 
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To  specify  the  type  of  output  desired  for  a  given  set  of  data^  assign 
the  proper  value  to  the  input  integer  variable  MD.   MD  can  take 
on  values  1  to  3  with  the  following  significance  attached  to  each: 
MD  =  1  Tabular  Output  only 

MD  =  2  Graphic  Output  only 

MD  =  3  Both  Graphic  and  Tabular  Output. 

If  MD  is  set  equal  to  1,  ten  columns  of  calculated  values  will 
be  printed  out  for  all  of  the  desired  combinations  of  secondary 
effects.   The  values  shown  are:   the  value  of  the  subscript  (N) , 
the  value  of  the  integer  MT,  the  combination  designator,  moment, 
moment  ratio,  the  converged  values  of  F-,  and  F^ ,  the  value  of  the 
parameter,  and  the  number  of  iterations  required  for  convergence 
in  each  case.   In  addition  at  the  top  of  each  page  of  output  are 
listed  the  values  of  the  raw  data  exactly  as  it  is  read  in  from  the 
data  card.   These  values  only  serve  as  a  convenience  in  showing  what 
was  read  into  the  program  and  what  the  individual  values  were  prior 
to  the  parameter  loading.   Other  values  that  are  printed  include  the 
exact  moment,  the  parameter  that  was  investigated,  the  elapsed  time 
in  seconds  since  the  last  stored  parameter  value  with  its  associated 
values  were  calculated,  and  the  critical  loads  for  each  parameter 
except  the  leg  ratio.   Actually  the  output  sheet  is  very  straightforward 
if  one  follows  the  column  labels. 

Some  of  the  above  mentioned  terms  deserve  a  little  more  explanation. 
Combination  Designator:   The  four  character  designator  supplied 
to  the  main  program  by  the  subroutine  Factor.   The  four  characters  are 
either  zero  or  one  and  they  represent  the  value  of  a  particular  factor 
for  the  given  combination  of  secondary  effects.   For  instance,  if  the 
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designator  is  1001,  the  factor  that  is  multiplied  by  the  term  represent- 
ing the  secondary  effect,  arc  and  chord  difference,  is  1.   The  factor 
that  is  multiplied  by  the  secondary  effect  term  for  shear  is  0.   The 
factor  multiplied  by  the  term  representing  axial  deformation  is  0  and 
the  factor  multiplying  the  beam  column  term  is  1.   The  order  of  the 
factors  always  appear  in  the  same  order.  Arc  and  chord,  shear,  axial 
deformation,  beam  column  action  appear  in  that  order  from  left  to  right. 

Exact  Moment:   When  the  term  exact  is  used  referring  to  a  calculated 
quantity  in  the  system  of  equations  it  means  that  the  calculation  was 
made  with  all  four  of  the  secondary  effects  allowed  and  none  of  the 
effects  suppressed;  for  such  a  calculation,  the  combination  designator 
is  1111. 

Moment  Ratio:   Ratio  of  moment  calculated  with  given  combination 
designator  to  exact  moment. 

If  MD  is  set  equal  to  2,  only  graphic  output  results.   For  graphic 
output  certain  limitations  exist: 

a.  The  ordinate  is  always  the  moment  ratio. 

b.  The  abscissa  is  always  the  parameter  under  investigation. 

c.  MA  or  the  number  of  stored  values  for  each  subscripted 
variable  should  not  exceed  30. 

One  graph  is  obtained  for  each  set  of  data  investigated.   On  each 
graph  six  sets  of  values  are  plotted.   These  six  sets  correspond  to  the 
first  six  values  of  MT.   See  Subroutine  Factor  in  Section  3.8  for  the 
combinations  of  secondary  effects  that  correspond  to  the  six  values. 
The  combinations  included  in  the  six  are  with  all  effects  allowed  (exact 
solution),  all  effects  suppressed,  four  combination  where  only  1  of 
the  four  effects  is  allowed  and  the  other  three  are  suppressed.   There- 
fore, if  the  user  specifies  graphic  output  only,  the  input  integer 
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variables  MB  and  MC  should  equal  to  1  and  6  respectively.   Calculations 
for  all  other  combinations  would  be  redundant  since  there  is  no 
tabular  output  and  the  first  six  combinations  are  all  that  will  be 
represented  on  the  graph.   If  on  the  other  hand,  tabular  output  is 
selected  or  the  user  desires  both  types  of  output  on  the  same  run 
then  the  above  restriction  is  not  applicable. 

Before  one  employs  the  graph  plot  form  of  output  it  is  necessary 
to  have  some  knowledge  of  the  values  that  are  to  be  plotted.   The 
data  card  contains  four  of  the  arguments  to  the  DRAW  subroutine. 
For  a  meaningful  graphical  output  a  reasonable  scale  has  to  be 
selected  and  other  parameters  have  to  be  chosen.   The  recommendation 
therefore  is  to  use  the  tabular  output  first  and  from  it  plan  the 
graphical  output  so  as  to  describe  the  results  most  usefully.   In 
addition  the  user  should  be  familiar  with  the  information  relating  to 
the  DRAW  subroutine  contained  in  the  next  section  and  in  the  subroutine 
system  library  write-up  reproduced  in  Appendix  A-4. 
3-8  Subroutines 

There  are  three  subroutines  used  in  the  program.   The  first  two 
subroutines,  designated  ALPHA  and  FACTOR  were  programmed  by  the  author 
and  the  third,  designated  DRAW,  is  from  the  subroutine  library  at 
the  computer  facility  of  the  Naval  Postgraduate  School,  Monterey, 
California.   The  complete  library  write-up  for  DRAW  and  its  use  appears 
in  Appendix  A-4. 

ALPHA:   This  subroutine  is  used  by  the  main  program  to  compute 
the  values  of  the  alpha  coefficients  as  they  are  shown  in  Table  II, 
page  13  .   It  is  called  twice  during  each  iteration  and  supplies  the 
program  with  values  of  the  C\S   or^jcorresponding  to  the  information 
that  is  given  to  it,  name  ly ,  &/ ,  /)  ,  and  EP1  or  7/S ,  f\    ,    and  EP1, 


n   ,  and  EP1  orW  ,r\ 
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depending  on  whether  the  main  program  wants  the  values  of  o(  S 
or  C{  S  . 

FACTOR:   This  subroutine  is  used  by  the  main  program  to  set 
the  combination  of  secondary  effects  to  be  considered.   As  previously 
mentioned,  it  works  along  with  the  integer  variable  MT  which  takes 
on  values  1  through  16.   Each  of  the  values  designates  a  certain 
combination  of  the  secondary  effects.   Upon  receiving  this  number 
designation,  the  subroutine  supplies  the  program  with  four  factors 
each  of  which  has  the  value  of  either  1.0  or  0.0  depending  on  the 
combination.   These  factors  are  multiplied  by  the  quantities  that 
introduce  the  secondary  effects  into  the  system  of  equations.   If 
the  factor  is  1.0  that  particular  secondary  effect  is  included  in 
the  calculations;  if  it  is  0.0  the  secondary  effect  is  suppressed  in 
the  calculation.   Tables  IV  and  V  on  the  following  two  pages  give  more 
information  concerning  these  combinations  and  their  designations 
that  is  helpful  in  understanding  the  particular  system  for  varying 
the  combinations. 

DRAW:   This  subroutine  is  the  one  through  which  the  program 
affords  the  user  the  graphic  output  option.   The  complete  write-up 
of  this  subroutine  just  as  it  appears  in  the  Postgraduate  School's 
subroutine  library  is  reproduced  in  Appendix  A-4.   It  is  felt  that 
no  further  explanation  is  necessary  except  for  the  four  arguments 

a 

that  are  read  in  as  data  in  the  main  program.   Referring  to  the 
write-up  in  Appendix  A-4,  the  arguments  h.  through  k.  are  supplied 
to  the  calling  statements  by  the  last  four  values  on  the  second 
data  card.   This  allows  user  control  of  the  scales  desired  and  the 
axis  location. 
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TABLE  IV 

FACTOR  NAME  SECONDARY  EFFECT  THAT  IT  CONTROLS 

FAC1  Difference  in  Arc  and  Chord 

FAC2  Shear 

FAC3  Axial  Deformation 

FAC4  Beam  Column  Action 
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SECTION  IV 
RESULTS  OF  INVESTIGATION 

4-1  Errors  in  Paper  No.  63--APMW-21 

Prior  to  developing  the  computer  program  to  utilize  the  relation- 
ships of  Section  II,  certain  review  and  checking  was  undertaken  by 
the  author.   This  review  was  undertaken  principally  to  gain  a  working 
familiarity  with  the  equations  and  nomenclature.   During  the  course 
of  this  review  and  later  in  developing  the  program  a  few  errors 
were  detected  in  the  original  paper. 

Most  of  these  mistakes  were  algebraic  in  nature.   One  was  in 
the  area  of  the  critical  buckling  loads  for  beam  column  considerations, 
The  mistakes  are: 

1.  Equation  (8),  the  factor  "2.0"  in  the  second  term  of  the 
numerator  did  not  appear  in  the  same  equation  in  Brock's  paper. 

There  were  three  mistakes  in  the  formulas  for  the  o(   coefficients 
of  Table  II. 

2.  &  3.   In  the  formula  for  Ok  ,  theG($   in  the  first  two  terms 
of  the  series  should  be  o(,5.   In  the  reference  paper  they  are  o(f  $  . 

4.  The  formula  for  of.±  ,    the  last  term  should  be  ffcLas 
opposed  to  Voj  "  -  the  paper.  This  was  apparently  a  typographical 
error. 

5.  The  last  incorrect  item  found  was  in  the  statement  relating 
to  the  critical  buckling  loads.   The  paper  indicates  that  2^/-  '* /A. 
defines  the  critical  buckling  load  for  the  ell-shaped  pipe  structure. 
This  is  true  only  for  the  equal  leg  case  but  is  not  true  in  general. 

All  of  the  results  reported  in  this  paper  reflect  the  corrections 
as  stated  above. 
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4-2  Numerical  Calculations 

With  the  use  of  the  digital  computer  program  described  in 
Section  III,  calculations  were  made  for  several  sizes  of  pipe 
and  for  various  values  of  temperature  increase,  as  well  as  for 
different  values  of  the  other  inputs.   This  was  done  not  only  to 
test  the  integrity  of  the  computer  program,  but  also  to  obtain 
sufficient  data  to  permit  a  search  for  trends  and  useful  modes 
of  presentation.   It  turns  out  that  the  leg  ratio,  a/b,  is  the 
most  significant  parameter,  as  might  have  been  expected,  and 
presentations  are  confined  to  displaying  the  moment  ratio  as  a 
function  of  leg  ratio  for  various  combinations  of  secondary 
effects  all  curves  shown  being  for  the  following  specific  case. 


Length  of  leg  a 
Length  of  leg  b 
pipe  inside  diameter 
pipe  outside  diameter 
Temperature  difference 

pipe  internal  pressure 

pipe  modulus  of  Elasticity 

pipe  coefficient  of  linear 
expansion 


varied 

360  inches 

10.020  inches 

10.750  inches 

530  degrees  fahrenheit 
(70-600  degrees  F) 


400  psi 

27.4  x  106  psi 

7.23  x  10"6  in/in  deg  F 


4-3  Final  Results  and  Conclusions 

The  results  are  shown  in  the  form  of  eight  curve  sheets, 
Figures  6  through  13,  which  should  be  considered  in  pairs,  the  first 
of  each  pair  showing  the  effect  of  suppressing  (i.e.  failing  to 
consider)  only  a  single  secondary  effect  and  the  second  showing  the 
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effect  of  suppressing  the  other  secondary  effects  and  considering 
only  this  one.   For  comparison  purposes,  each  curve  sheet  also 
shows  the  effect  of  neglecting  all  four  secondary  effects. 

These  curves  are  based  only  upon  what  has  been  called  the 
"moment  ratio1',  i.e.,  the  ratio  of  the  moment  at  the  intersection 
of  the  two  legs  as  calculated  by  an  inferior  procedure  to  that 
calculated  by  the  best  available  procedure,  namely  that  of  including 
all  four  secondary  effects.   It  is  clear  that  substantially  the  same 
sort  of  result  would  have  been  obtained  by  considering  a  similar 
moment  ratio  for  either  of  the  two  built-in  ends  or  a  similar  ratio 
of  reaction  force. 

These  curves  may  be  regarded  as  speaking  for  themselves.   The 
most  significant  conclusion  that  may  be  made  is  that  neglecting 
secondary  effects,  indeed   can  have  a  profound  influence  on  the 
accuracy  of  piping  flexibility  calculations,  leading  to  gross  errors 
in  the  evaluation  of  stresses  in  cases  of  leg  ratios  which,  though 
extreme,  are  likely  to  be  encountered  in  piping  practice. 

Since  the  secondary  effects  are  quite  subtle  to  evaluate,  this 
conclusion  could  not  have  been  forthcoming  from  a  less  sophisticated 
and  demanding  approach.   It  is  believed  to  be  an  entirely  new  con<- 
clusion  and  it  is  certainly  a  significant  one. 
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APPENDIX  A 
COMPUTER  PROGRAM  NOMENCLATURE 


A-l  MAIN  PROGRAM 

NOTE:  Where  possible  the  variable  names  of  the  computer  program 
are  identified  by  the  notation  of  Tables  I  and  II.  Also,  in  some 
cases  the  same  physical  quantity  may  be  represented  in  the  program 
by  more  than  one  variable  name.   In  such  cases  the  extreme  right 
hand  column  will  indicate  the  general  area  of  the  program  where  the 
physical  quantity  takes  on  the  indicated  variable  name. 


VARIABLE  NAME 

A 

AA(N) 

A1-A10 

AB1-AB10 

ALIT(N) 

ALPHA 

AREA 

B 

BB(N) 

Bl 

B2 

BET1 

BET2 

COEFI 

DELT1-DELT4 

DET 

DIN 

DOUT 


TABLE  I  &  II  NOTATION 
a 
a 

oc,  —  o</0 

e 
Subroutine  ALPHA 

A 
b 
b 


COMMENT 
Computation  Section 
Loading  Section 


Loading  Section 


Computation  Section 
Loading  Section 


feff>H)b-  J 
[a'fa,  +  ^)]/(aaj -t-  bol3) 
Cbl  ft +ij^]/(CL«i+  bW) 


Material  coefficient  of  thermal  expansion 

Size  of  6par  for  each  parameter  1-4  respectively 


pipe  inside  diameter 
pipe  outside  diameter 


Computation  Section 
Computation  Section 
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COMPUTER  PROGRAM  NOMENCLATURE  (cont.) 


VARIABLE  NAME 

DTEMP 

E 

ELIT 

EP1,EP2 

ETA 

ETAB 

EXMOM(N) 

EXSCAL 
FAC1-FAC4 

FACTOR 
F1,F2 

F1CRIT,F2CRIT 
F1N,F2N 

FTE11-FTE162 
G 

H 

I 

IA 

LABEL 

ITITLE 


TABLE   I  &   II  NOTATION 


COMMENT 


Difference  of  Temperature 


Computation  Section 


*»* 


0 
5 


Exact  moment  for  a  particular  parameter  variable,  i.e. 

with  all  secondary  effects  allowed. 
Applicable  to  graphic  output:   X  scale  in  units  per  inch 
Four  factors  used  for  either  suppressing  or  allowing 
their  corresponding  secondary  effects. 
Subroutine  FACTOR 
Assumed  values  of  the  forces  F-|>F9  at  the  beginning 
of  an  iteration. 

Euler  Buckling  load  for  F  and  F  . 
Newly  calculated  values  of  F-  and  F  from  a  particular 
iteration. 

See  section  on  input  Data 

G  Computation  Section 

coefficients  defined  in  equations  (16)- (19) 
H 
Counter  use  to  count  the  number  of  iterations  required 

for  convergence. 
Factor  to  aid  in  determining  improved  values  of  F's 
Argument  for  Subroutine  DRAW 
Argument  for  Subroutine  DRAW 
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COMPUTER  PROGRAM  NOMENCLATURE  (cont.) 


VARIABLE  NAME 


TABLE  I  &  II  NOTATION 


COMMENT 


IXUP 

IYRIGH 

J 

KHAT 

KKK 

KPAR1 

L 

LABEL 

LL1 

MA 

MB 

MC 

MD 

MIT 


MT 

N 

NN 

P 

PAR2 


Argument  for  Subroutine  DRAW 
Argument  for  Subroutine  DRAW 
Counter  used  to  count  the  number  of  subincrements  between 

the  larger  main  increment. 
Secondary  Effect  combination  designator 
Input  variable  that  assigns  value  to  IA 

Input  variable  that  designates  parameter  to  be  investigated 
Input  variable  that  aids  in  determining  ^fcPAR. 
Argument  for  Subroutine  DRAW 

Input  variable  that  sets  size  of  subincrements,  opar. 
Input  variable  that  designates  final  value  of  subscript, N. 
Initial  value  of  MT 
Final  value  of  MT 

Input  variable  that  designates  type  of  output 
Input  variable  that  specifies  the  maximum  number  of  iterations 

that  will  be  performed  before  terminating  operation  on  a 

particular  set  of  data. 
Integer  variable  that  specifies  the  combination  of  secondary 

effects  to  be  considered. 
Subscript,  used  to  identify  calculated  values  of  desired 

quantities  with  the  corresponding  parameter  value. 
Initial  value  of  the  subscript  N 

p  Computation  Section 

A  designation  for  the  second  parameter,  (Radius  of 

Gyration) /(Total  Length),  in  the  computational  part 

of  the  program. 
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COMPUTER  PROGRAM  NOMENCLATURE  (cont.) 


VARIABLE  NAME 

PIE 

PIESQ 

POIS 

Q 

R 

RAGY 

RHO 

S 

SFAC 

W 

WBAR 

XDIN(N) 

XDOUT (N) 

XLA 

XLAB 

XMU1 

XMU2 

XN 

XP(N) 

Y  SCALE 

XI 

ZIN 

ZLEGR(N) 

ZLEN 

ZMOM 


TABLE  I  &  II  NOTATION 

7T 


COMMENT 


■J 


Radius  of  Gyration 


Computation  Section 
Input  Data  Section 
Input  Data  Section 


Input  Data  Section 


Loading  Section 
Loading  Section 


f 

b 

? 
If 

w 

pipe  inside  diameter 
pipe  outside  diameter 

fa 

Moment  in  Computation  Section 

p  Loading  Section 

Applicable  to  graphic  output:   Y  scale  in  units  per  inch 

I 
pipe  inside  diameter         Input  Data  Section 
First  parameter,  a/b 
Total  Length,  a  +  b 
moment 
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COMPUTER  PROGRAM  NOMENCLATURE  (cont.) 


VARIABLE  NAME     TABLE  I  &  II  NOTATION 


COMMENT 


ZOOM  Initial  thermal  strain 

ZOUT  pipe  outside  diameter 

ZPAR2(N)  Second  parameter,  (Radius  of  Gyration) /(Total  Length) 

ZPAR3(N)     Third  parameter,  Pressure /Modulus  of  Elasticity 

ZRAT1(N)-ZRAT16(N)   Ratios  of  the  Moments  of  the  sixteen  combinations 

divided  by  the  exact  moment  for  that  combination. 
ZRATIO       Same  as  the  ZRATQ'S  except  this  value  is  not  stored. 


A- 2   SUBROUTINE  ALPHA 

VARIABLE  NAME  TABLE   I   &  II  NOTATION 


c 

?\  or  r\ 

EP1 

6/ 

T1-T10 

<*/  "  <X,o   °r     °<,-i 

TX 

iwl   or  /©/ 

Y 

l-isfor  (-&)* 

Z 

(Vpor  flj)*- 

ct 


/o 


COMMENT 


A- 3  SUBROUTINE  FACTOR 

VARIABLE  NAME     TABLE  I  &  II  NOTATION 


COMMENT 


JT  Integer  variable,  MT 

LLAT         Secondary  effect  combination  designator 
R1-R4        Four  factors  used  for  either  suppressing  or  allowing 
their  corresponding  secondary  effects. 
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A- 4   SUBROUTINE  DRAW-DESCRIPTION  FROM  SUBROUTINE  LIBRARY 

A.  IDENTIFICATION: 

TITLE:  General  Graph  Output  Subroutine 

SUBROUTINE  NAME:   DRAW (J7-NPS- DRAW  for  CDC  1604  revised  for  IBM  360) 

DATE:  July,  1967 

PROGRAMMER:       Revised  by  Patricia  Johnson 

B.  PURPOSE: 

This  subroutine  generates  data  for  plotting  graphs  on  a  CalComp 
plotter.   Provision  is  made  for  curve  drawing  and  point  plotting, 
automatic  scaling,  graph  titling  and  axis  annotating.   An  attempt 
was  made  to  provide  a  considerable  amount  of  flexibility,  at  the 
expense,  necessarily,  of  a  relatively  large  number  of  arguments 
and  a  rather  high  memory  requirement. 

C.  USAGE: 

1.  Definitions: 

In  what  follows  the  word  "graph"  will  be  taken  to  mean  one 
piece  or  frame  of  graph  paper  on  which  there  may  be  plotted 
one  or  more  curves  and/or  sets  of  points.  A  "curve"  will 
mean  a  continuous  line  generated  by  the  sequence  of  straight 
lines  joining  successive  points  of  the  set  defining  the 
curve.   A  "point  plot"  will  describe  the  representation  of 
a  succession  of  points  by  means  of  symbols  (such  as  a  cross) 
on  the  graph.   The  points  are  not  connected  in  a  point  plot. 

2.  Calling  Arguments: 

All  necessary  information  is  transferred  to  DRAW  through  the 
calling  arguments.   The  call  statement  is:   CALL  DRAW  (NUMPTS, 
X , Y ,M0DCUR , ITYPE , LABEL , ITITLE ,EXSCAL , YSCALE , IXUP , IYRIGH , 
MODXAX ,M0DYAX , IWIDE , IHIGH , IGRID ,LAST) 
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It  is  important  to  realize  that  one  and  only  one  curve  or 
set  of  points  is  plotted  each  time  DRAW  is  called.   However, 
it  is  possible  to  call  DRAW  repeatedly  if  several  curves 
and/or  sets  of  points  are  wanted  on  one  graph. 
The  calling  arguments  are  as  follows: 

a.  NUMPTS:    The  number  of  points  defining  a  curve  (21 NUMPTS S  900) , 

of  the  number  of  points  to  be  point  plotted  (2S  NUMPTS130) 

b.  X:         The  array  of  X-ordinates,  X,  must  be  dimensioned  at 

least  equal  to  NUMPTS  in  the  calling  program. 

c.  Y:         The  array  of  Y-ordinates,  Y,  must  be  dimensioned  at 

least  equal  to  NUMPTS  in  the  calling  program. 

Controls  the  number  of  curves,  and/or  sets  of  points 

on  one  graph: 

This  is  the  only  curve,  or  set  of  points,  to  be 

plotted  on  this  graph. 
=  1         This  is  the  first  of  two  or  more  curves,  and/or  sets 

of  points,  to  be  plotted  on  this  graph. 
"  2         This  is  an  intermediate  curve,  or  set  of  points. 
=  3         This  is  the  last  curve,  or  set  of  points,  for  this  graph. 
e.   ITYPE:     Controls  the  type  of  plot  (i.e.,  curve  or  point  plot): 
=  0         This  set  of  points  is  to  be  represented  by  a  curve . 
=  1         These  points  are  to  be  plotted  with  a  cross  (x) . 
=  2         These  points  are  to  be  plotted  with  a  plus  (+) . 
=  3         These  points  are  to  be  plotted  with  a  square  (□) . 
=  4         These  points  are  to  be  plotted  with  a  diamond  (O). 
=  5         These  points  are  to  be  plotted  with  a  triangle  (A). 


d.  MODCUR: 


=  0 


59 


f.   LABEL:     This  is  a  Hollerith  curve  or  point  identifier. 

If  a  curve  is  being  drawn,  LABEL  must  have  4  char- 
acters (including  any  blanks),  and  these  will  be 
reproduced  beside  the  end  of  the  curve.   This 
argument  can  be  set  in  the  calling  program  by  a 
declaration  statement,  e.g., 
REAL  LABEL/4H1234/ 
or        REAL  LABEL/ f 1234' / 
or        REAL  LABEL/4H     /    (The  latter  must  be 

used  when  no  label  is 
wanted.) 
If  a  set  of  points  is  being  plotted,  LABEL  is  an  8- 
character  identifier.   The  first  4  characters  will  be 
reproduced  beside  the  first  point,  and  the  last  four 
characters  will  be  reproduced  alongside  the  last 
point.   This  argument  can  be  set  by  a  declaration 
statement,  e.g., 

REAL*8   LABEL/8HFRSTLAST/ 
or        REAL  LABEL(2) /4HFRST,4HLAST/ 
or        REAL*8  LABEL/8H  /   (The  latter 

must  be  used 
when  no  label 
is  wanted) . 
The  above  arguments,  a.  through  f.  (and  q.),  have 
meaning  every  time  DRAW  is  called.   On  the  other 
hand,  the  remaining  arguments,  g.  through  p.,  have  no 
meaning  except  when  MODCUR  =  0  or  1. 
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g.   ITITLE:    An  array  of  96  EBCDIC  characters,  the  first  48  of 
which  will  form  the  first  title  line,  and  the  last 
48  the  second.   The  array  must  be  dimensioned  in 
the  call ins  program,  must  contain  the  user's  job 
identification,  and  must  have  unwanted  characters 
set  to  blank.   This  may  be  done  with  a  declaration 
statement  such  as: 

REAL*8  ITITLE (12) /'96EBCDIC  character'/ 

h.   EXSCAL:    X- scale  in  units  per  inch.   EXSCAL  will  always  be 

rounded  off  to  one  figure  significance.   If  EXSCAL  =  0, 
the  X-scale  will  be  computed  by  DRAW.   This  is  called 
auto-scale. 

i.  YSCALE :    Y-scale  in  units  per  inch,  with  properties  corresponding 
to  those  of  EXSCAL. 

j.   IXUP:      Distance,  in  inches,  of  the  X-axis  from  the  bottom  of 
the  graph  (O^IXUP  ^IHIGH)  .   This  argument  will  be 
ignored  unless  MODXAX  =  2,  see  below. 

k.  IYRIGH:  Distance,  in  inches,  of  the  Y-axis  from  the  left  of 
the  graph  (0*IYRIGH£IWIDE) .  This  will  be  ignored 
unless  MODYAX  =  2,  see  below. 

1.   MODXAX:    Determines  the  mode  of  the  X-axis  location: 

=  0         The  X-axis  will  be  located  automatically  by  DRAW, 

with  the  origin  of  Y  on  the  graph. 
=  1         The  X-axis  will  be  automatically  located  by  DRAW, 

with  the  origin  of  Y  moved  (in  one's  imagination)  an 
integer  number  of  inches  above  or  below  the  graph,  if 
this  is  appropriate.   This  option  can  be  used  only  if 
the  Y-scaling  is  automatic  (YSCALE  =  0). 
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=  2 
m.   MODYAX: 

n.   IWIDE: 

o.   IHIGH: 

p.   IGRID: 


q.   LAST: 

=  0 
=  1 
=  2 


=  3 


The  X-axis  location  will  be  as  specified  by  IXUP. 

Determines  the  mode  of  Y-axis  location  in  the  same 

way  as  MODXAX,  above,  governs  the  X-axis  location. 

Width  of  graph  in  inches  (1:£  IWIDES  9)  .   If  IWIDE  is 

out  of  this  range,  a  value  of  8  will  be  assumed. 

Height  of  graph  in  inches  (IS  IHIGH<  15)  .   If  IHIGH 

is  out  of  range,  a  value  of  8  will  be  assumed. 

If  IGRID  =1,  a  1"  x  1"  grid  will  be  superimposed  on 

the  graph.   This  is  useful  only  if  plain  paper  is 

used  on  the  CalComp  Plotter. 

Indicates  to  the  calling  program  whether  the  previous 

plot  was  completed  successfully.   The  codes  are: 

Last  plot  was  completed  successfully 

Last  plot  was  not  completed  successfully. 

Last  plot  was  not  completed  successfully  and  no  further 

graph  output  will  be  attempted  until  DRAW  is  next 

entered  with  MODCUR  =  1  or  0 

An  attempt  was  made  to  enter  DRAW  with  MODCUR  +   1  or  0 

while  the  error  lockout  was  set.   This  argument  must 

always  be  a  variable  name  and  never  a  number  in  the 

call  statement. 
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3.   Notes  and  Comments : 

a.  The  graph  scales  and,  if  MODXAX  =  1  and/or  MODYAX  =  1,  the 
amounts  of  origin  offset  are  always  output  as  part  of  the 
graph  title. 

b.  There  are  internal  checks  of  the  input  to  DRAW  to  prevent 
incorrect  use.   If  an  input  error  is  detected,  an  attempt 
will  be  made,  where  possible,  to  complete  the  plot.   If 
an  argument  is  "corrected"  in  this  process,  the  user  will 

be  so  informed  on  the  standard  output.   If  it  is  not  possible 
to  complete  the  plot,  the  user  will  be  informed  of  the 
reason  by  a  message  on  the  standard  output. 

c.  If  part  or  all  of  a  curve  would  fall  more  than  0.6"  laterally 
beyond  the  ends  of  the  X-axis,  or  0.7"  vertically  beyond 

the  ends  of  the  Y-axis,  the  X  and/or  Y  ordinates  will  be 
limited  so  that  the  curve  will  typically  become  a  line  along 
part  or  all  of  the  boundary  of  the  graph  as  here  defined. 

d.  If  one  or  more  points  of  a  point  plot  would  fall  outside  the 
graph  area,  the  plot  of  the  point,  or  points  will  be  inhibited, 
The  number  of  such  points  will  be  reported  to  the  user  on  the 
standard  output. 

e.  It  should  be  pointed  out  that  the  X  and  Y  scaling  and  axis 
locating  processes  are  entirely  independent,  so  that,  for 
example,  X  might  be  auto-scaled,  while  the  Y-scale  is 
specified.   At  the  same  time  the  X-axis  might  be  located 
automatically,  while  the  Y-axis  location  is  specified. 

f.  It  must  be  remembered  that  the  scales  and  axis  locations  of 
a  multi-plot  graph  are  set  when  DRAW  is  called  for  the  first 
time  (with  MODCUR  =  1).   Thus  the  user  must  attempt,  at  the 
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time  to  achieve  scaling  and  axis  location  which  will  be 
appropriate  to  all  the  plots  he  intends  to  make  on  the 
one  graph.   Particularly  if  the  automatic  features  of 
DRAW  are  selected,  foresight  will  be  demanded  of  the  user 
in  this  respect. 

4.  Auto-Scale  Properties: 

The  scale  factor  is  chosen  from  amongst  the  values  1,  2,  or  5 
units  per  inch,  or  some  power  of  10  times  one  of  those  values. 
A  curve,  or  set  of  points  which  is  plotted  with  auto-scale  will 
normally  lie  entirely  within  the  graph  area  is  defined  in  3.d., 
above.   The  only  exception  may  occur  if  an  axis  is  placed,  by 
the  user,  along  one  edge  of  the  graph  (e.g.,  IXUP  =  0,  MODXAX  =  2) 
In  such  a  case,  points  "outside"  the  axis  are  not  considered  in 
the  selection  of  a  scale  factor  (e.g.,  negative  X.  do  not  affect 
the  choice  of  scale  when  IXUP  =0).   If  automatic  axis  location 
as  well  as  auto-scale  is  selected,  the  plot,  if  it  does  not  fill 
the  graph  area,  will  be  placed  as  far  as  possible  towards  the 
bottom- left  of  the  graph  area  consistent  with  the  fact  that  the 
axes  can  be  set  in  increments  of  1"  only. 
1st  revision  -  31  July  1967. 
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APPENDIX  B 
COMPUTER  PROGRAM  FLOW  CHART 


START 


DIMENSION 


I 


(  READ  FIRST  SET  OF   DATA  ) 

4 

AsE-0?\ 


YES 


GO   TO  END 


[n< 


LOAD  PARAMETER 


INITIALIZE 


1 


SUBROUTINE  FACTOR 

Select  factors 
for  secondary  ef- 
fect combination 


PRELIMINARY  CALCULATIONS 


ASSUME  INITIAL  VALUES  OF  F'S 


COMPUTE  W  AND  W 


SUBROUTINE  ALPHA 
Compute  cVTs  orof's 


COMPUTE  NEW  VALUES   OF  F'S 


I 


TEST  OLD  F'S  AGAINST  NEW  F'S 


I 


IMPROVE 
VALUES 
OF  F'S 


IS  CONVER- 
-22— (    GENCE  CRITERIA 
SATISFIED  ? 


Y6S 


® 


© 


<S)   © 
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COMPUTE  MOMENT  ETC,  USING 
CONVERGED  VALUES  OF  F'S 


CHANGE 
COMBINATION 


ASSIGN  CONVERGED  VALUES  OF 
F'S  TO  PROPER  FTE 


NO 


IS   TABULAR 
OUTPUT  SPECIFIED   ? 


7y*\ 


C  PRINT  RESULTS  ) 


HAVE  ALL  COMBI-  \ 
NATIONS  OF  SECON- 
D1RY  EFFECTS  BEEN 
INVESTIGATED  ? 


YES 


HAS  ENTIRE  RANGE 
OF  PARAMETER  BEEN 
INVESTIGATED  ? 

lyes 


Mo 


©   © 


INCREMENT 
PARAMETER 


IS  GRAPHICAL  OUT- 
PUT SPECIFIED 


Yes 

(  CALL  DRAW) 


UT-\ 


END 


66 


APPENDIX  C 
DATA  CARD  DESCRIPTION 

C-l  Data  Card  Placement 

For  one  set  of  data,  six  data  cards  are  required  to  give  the 
program  the  necessary  input.   The  first  card  gives  the  physical  data 
such  as  leg  length,  pipe  size  and  material,  and  other  quantities  that 
describe  the  physical  situation.   The  second  card  supplies  the  control 
variables  to  the  program.   These  control  variables  indicate  to  the 
program  what  parameter  is  to  be  investigated  and  its  range. 

The  last  four  cards  supply  the  initial  values  given  to  the 
FTE's  (see  page  26).   Each  of  the  last  four  cards  contain  eight  of 
these  values.   Unless  a  particular  computer  run  starts  near  the 
critical  part  of  the  parameter  range  (where  convergence  is  difficult) 
the  values  on  the  last  four  cards  are  usually  zero.   If  such  is  the 
case,  blank  cards  may  be  used;  however,  when  values  other  than  zero 
are  supplied  as  input,  the  format  shown  in  Appendix  C-2  must  be  used. 

It  must  be  remembered  that  six  data  cards  must  always  be  used 
for  a  particular  set  of  input  whether  or  not  the  FTE's  are  zero. 

Any  number  of  sets  of  six  data  cards  can  be  supplied  to  the 
program  for  any  one  run.   Consideration  should  be  given  to  time 
requirements,  however.   After  the  final  set  of  six  data  cards,  one 
blank  card  follows  and  then,  after  that,  the  last  control  card  (see 
Appendix  D) .   The  blank  card  is  used  to  terminate  the  run  by  telling 
the  program  that  all  data  has  been  considered. 

The  user  should  be  familiar  with  fixed  and  floating  point  variables 

3 
and  the  common  Fortran  field  specifications. 


3  Daniel  D.  McCracken,  A  Guide  to  Fortran  IV  Programming  (New  York: 
John  Wiley  and  Sons,  Inc.,  1965),  pp.  3-5,  13,  84. 
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Appendix  C-2  shows  a  typical  input  (all  six  data  cards).   The 
extreme  left  hand  column  indicates  the  data  card  number  (this  order 
must  always  be  followed).   For  each  data  card  there  are  two  rows 
shown.   The  top  row  indicates  the  variable  name  whose  value  is 
shown  directly  below  in  the  next  row.   This  second  row  is  the 
actual  value  of  the  variable  named  above  it  and  appears  just  as  the 
data  card  itself  does.   The  top  row  of  the  two  in  each  case  serves 
only  to  aid  in  identification  of  the  variable. 
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C-2     Sample  Data  Card  Set 
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Blank  Caud 

onbset    I    &5 

OF  *4 

SIX 
DArA 
CA&bS 


//PIPEELL„^Jo6^'0\4t*\'JM.TMl&')MiGLE\/£L=L  \ 


APPENDIX  D 
CONTROL  CARD  ARRANGEMENT 
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APPENDIX  E 
PROGRAM  LISTING 


E-l  Main  Program 


X<       <N, 

h-a     uc     3r        <a  ■-»— <x            xr  <oox  >  rvrv  ».>o 

<■—<      5"  <J      2       ooo  wt-OK                X  acLLh-  *  -4-00— «—< 

DO.          a      O.      •—  •—  C  icKi^is.              cc  o>  _j  lulucjll 

o         luC     <     oocyuj  —or ^  »  czy  <  t— t— r-»h- 

UJO       Ik         I         >  U  Q-  CM       —  <*•—  t—  «_;  LLOLLULL 

U_'      k-oo      f*"      _l5  <  X    ►~-C— »  Q.       o  I/!  *  •■»—    ► 

u.a.        c     so     <xx  •—~>C'i-'0                                    u_  x  <-i»~iu.r\i 

C<r     U.Q.              2T2Too   •  -~r>i—  ,-4r-i         cos^  luZ_j  a.  -d-co  »lp 

X      O           •     <ir    i  r-  o— ><     — <         or-  xilu.  ••  lulupj  — 

5To~            _J       C                 _J>T  —i^a  —                 CC  L- LLOO  fV  h- I—.-1LU 

LU    I         V.<       Z       11    »       P  — I       N-<«-            — <LL                              a  CL  LLLL  — *k- 

(— _jO>                  X  on  LL  — >  — ■    *■— I O                  at  OK  LU  •■    -LL1  Li. 

oo_i     C<:     ex     i— >— c(~j  w-fr,— i-H             -5  oti/)<— "  *  r\r\it—  * 

>-lu     ac«t      iu        t/i       •■  a.iO<t-             K  u.<  —i  mMi-j 

oo         cc         a      olji/'Z  ck— 'cf<         coo  niO  _j  lull.'  »-tri 

z        o~     <     (— x*-1^  u.c.-^f^tx        ex                       a:  _j  »— I—  -4H 

LU<T          »LU       CL             k-OOX  _JNsl          »-M             -5Cf  OOZLL  •>  LLU_»-iLL' 

x        ilh—           5"     >  Nsi  •— » —  •■           k-  j^Ni  ^r  •>  -r-it— 

t— u.         <     oo     <_i_j  •»  •.•^.tro  —             oo  <r ex  ^  .—t— *luu. 

C      ZK      <      orO<TO  ~C-i-r-to                      w  3     1-  :*:  rof^k-  ► 

LU             XoO                  COZU  C— •  <T. — I. — >                  LU       LL  OCCS  ••  LULLLLt\J 

>o~     O         on     GKO  — '.-icl     — <             a      lu  lullz  t—  h-t-   •">!■ 

_i«->     -sCorao  h     n-              «i_     q             |_  _  LLu_rv— i 

OlO         uj     <i     owv  »  —.  »-o—         h- co  —  Xoo  S"  •■  -OLU 

O0>-        CLk-        LU                   k-OO  <— r\J  — — »LO.                •                »  l/;a>-     •  •■  <^jfV  . — 1»— 

_J       CJ>-h       Q_     •CCLU»-«  LU  5"DtOt—  ■"-•             XQ       LU  X5T       Z  O  <MOLULL 

X<     «OOZ        Q.  OO  LL1  k- _j  _J  C<T*- '<Tt—                 »LU          •■  jClL'C  7"  LULUk-    » 

t—zocoox     <ro>—  <=r  • — •■ — i  :fcl— >cc<-         to.     o  ooX)h->  •>  t— t—  ll--< 

<LULU                  "D^IJ  XNJ       NQt                   <         •■  C_J»—  O  LULL    »>*■ 

C       ZTULUJ    •O02T0LC*-'  LU    ••—    »-rvJ               ..ICL  CK<<  2T  »   ►Hr- 

Lu>ac:crx<'-,<^<x  •  •.— .>t-»  ►       oooo     s:  s:z>c!  -  r-nou 

ZhUivu-uiicaiik  — o>— O—         ail      uj             lu     •—  CO  oovOf-"!— 

0«—  UQ.       Z>OUO-J    •  Of-<^0            -J— i       ►—  t/>Q£UJI—  2:  LULULULL 

•—  _i         u_nt:_iLj     kllx  — < --•  ci ^ ^-i         »-i         a  -  orxon  •>  h-»—  \—  - 

00w-LU>OO<y  on i/l  — «       KsJ       —i       \5"u.          •>  OXt— a  <  — LLLULLfNJ 

LUfrarcr.     u.7     >-czz  is!     —  •■—                o     t—  zu     >  —ix  f\)  ►  -  -cr 

cji_<     h-<oia>-i  i  wf-.-'Cr-^                   3  3     u_2  r-  ••  >—  r\jrM(\i^-t     m 

xDllZj     u-t-         c_  ceo  — ot~vt             oo     c-~  Olu«—  •—  mz  (\j--loctlu     m 

OOLUOOU-KLJ.— ■       _„'oO-3  ►C_J— <<^h                  i-tr-^isji^  >X  -2i  -LULULLh-       r>- 

—  ^ioo<5:ox_jLL<t—  x  ►<— *zl\-         •  ooo  ►  •  «*i—     a:  cr  •■  >ot— t— h-ll      ••            _j 

LL       5"t-       (— O-OU;    •  I   —    ►       N<V       \>    iZO                     «C  0_'  •LLLLLL    »       — <  "V 

i_     <t     a:  *-     a.     <ujk-  <i"»'  -c£.          "_)•«■«—<— i  Luaca.  rvi  »  cc   ••  •>  »-^     m             ct 

<<r     i/i<>LL<zzy — i  -^o<^~  fvi         rv<«->^slC         cz  •»-"  ll-h— i— iro     r-  i 

Cr:       XOaUJCJ       O       LL  a'r-i-j(-C    ►      X^-^ZO    »-f\l  LLLLLUQi  U"  cc  m  ^iu>  cr  — i         » 

0?L  3:OJiX        U."— •     »LU>  ■«•         —  <— 4— .       CO-—  <I~^oO    ••  C              LU  G«I  -LL  LJ  LU  LU  — O                   — - 

LJU*-1        CjLUK-LJH- OO        QQ  _i^-        a-^UINUJ        LU    *«J-              <l_lt—  tMQ.  1^- 1— t— t— H— U  .  O  ^ 

o:cl-«oo     h- a.     <i  •>-  «g<r  — kj     i— i vT _j _j >■  2. ac  •  u.f—  t—  lu  i<c  •— 'Lullllll   »(nj            ~» 

Q-LL       wOZOwODOC    •■  ^J<Q_    •■—  -tS  LU  I—  ►—•-•— >0  Z3<I       5l  — — '          •■■ ^-.-^o—                  <X 

200^t_jQ.<--      <:  ct      x— co     _j  ra  •—  •—  l-  o  — •  _!Qoo<r  o^  Lnroromrn.--.-*              < 

oo'j>a>->—  2: luz_z2oj  zz  "Oh- — Lu<h-_j— -oll  <     lux  •  r- X >-" Is- h- r^  1*- Q  I      CtOO  II  o 
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E-2   Subroutine  ALPHA 
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E-3   Subroutine  FACTOR 
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E-4     Sample  Output 
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APPENDIX  F 
SAMPLE  CALCULATION 


The  following  table  of  values  is  provided  to  aid  the  user  in 
checking  the  computer  calculations.   Two  combinations  of  secondary 
effects  are  listed.   The  first  is  the  "Exact  Value"  combination 
where  all  secondary  effects  are  included.   The  second  is  the  most 
erroneous  combination,  or  the  one  where  all  secondary  effects  are 
suppressed  in  the  calculation.   The  descriptive  data  is  that  listed 
on  page  42  of  Section  IV  and  also  exactly  the  same  as  the  sample 
data  card  in  Appendix  C,  with  a  leg  ratio  of  0.2. 


Quantity 
Assumed  Fl 
Assumed  F2 
W 

«*, 
<** 
oh 

*l 

<*S 

<*io 

WBAE 

a, 

of, 


All  Effects  Allowed     All  Effects  Suppressed 


65372.66553 
4623.212071 
.0054419338 
.5011360506 
.3340605064 
1.001817935 
1.002727150 
.4185170188 
.2678453147 
.0072796731 
.6695786688 
.4320902252 
.2859409726 
1.923736751 
2.324235399 
1.496163752 


74373.10291 

5192.148178 

0.0 

.5000000000 

.3333333135 

1.000000000 

1.000000000 

.4166666865 

.2666666508 

0.0 

.6666666865 

.4166666865 

.2666666508 

0.0 

.5000000000 

.3333333135 
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Quantity 

All  Effects  Allowed 

All  Effects  Suppressed 

*» 

3.878225194 

1.000000000 

5(4 

5.471216862 

1.000000000 

as 

8.613517662 

.4166666865 

*<. 

5.485230450 

.2666666508 

0(-J 

31.81243915 

0.0 

#e 

13.54446691 

.6666666865 

5ed 

8.624985666 

.4166666865 

O<io 

5.499880474 

.2666666508 

* 

1.817163214 

.6000000000 

fc 

205.4115482 

150.0000000 

4 

-.0001242864 

-.0000882817 

<k 

.0018232014 

.0013242253 

&. 

.0000294904 

.0000247189 

4» 

-.0001211686 

-.8828170315 

A 

.0132454070 

0.0 

A 

.0132186939 

0.0 

Fl  Converged 

65840.74973 

74972.17486 

F2  Converged 

4633.905400 

5208.593999 

Moment 

1071501.071 

1231122.149 
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APPENDIX  G 
SHORTENING  DUE  TO  DIFFERENCE  BETWEEN  ARC  AND  CHORD 

The  secondary  effect  produced  because  of  difference  between 
arc  and  chord  is  described  by  equation  (7)  of  Section  II  when 
higher  order  terms  are  neglected.   The  following  is  a  more  detailed 
development  of  that  relationship. 

In  order  to  have  a  physical  understanding  of  the  background  of 
the  derivation,  refer  to  Fig.  G-l.   Due  to  the  action  of  cross  shear, 
as  it  is  idealized  in  customary  theory,  the  horizontal  length,  dx, 
from  0  to  A, is  stretched  as  it  is  inclined  to  angle  0  >  so  that  it  goes 
from  0  to  B.   Then  the  action  of  bending  further  inclines  it  so  that 
it  takes  the  position  OC ,  but  there  is  no  change  in  length  during  this 
bending  process.   (Change  in  length  due  to  axial  force  is  accounted 
for  separately.)   The  increment  of  difference  between  arc  and  chord 
is  the  distance  dA  shown  in  the  figure. 

In  the  derivation,  it  is  first  assumed  that  0  is  small  so  that 
tan  ft  =o  approximately. 


dy 


•   • 


ds2  =    dx*  +-   dy2, 
dh  =   dx-  ds  cos(y') 

^  -'  If  ^    dy«XdX 
ds*  =  dx*  *  *W 

d$*  dx'0+**) 

ds  a    dx  Vl-tY* 
dh-  dx-dxJTT^CosCy') 

d*=  dxO-fi+FcosCvj')) 
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If  series  representation  of  cos(U)  and  V/"^)f   i-s  substituted 
with  higher  order  terms  neglected 


or 


(7)      A  =  i    <J  [(yf  -  *l]dX 
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